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We consider single-particle spectra of a symmetric narrow-band Anderson impurity model, where 
the host bandwidth D is small compared to the hybridization strength Aq. Simple 2nd order 
perturbation theory (2PT) in U is found to produce a rich spectral structure, that leads to rather 
good agreement with extant Lanczos results and offers a transparent picture of the underlying 
physics. It also leads naturally to two distinct regimes of spectral behaviour, AqZ/D ^ 1 and <^ 1 
(with Z the quasi-particle weight), whose existence and essential characteristics are discussed and 
shown to be independent of 2PT itself. The self-energy Ei (lj) is also examined beyond the confines 
of PT. It is argued that on frequency scales of order lu ~ \/ AqD, the self-energy in strong couphng 
is given precisely by the 2PT result, and we point out that the resultant poles in Ei((j;) connect 
continuously to that characteristic of the atomic limit. This in turn offers a natural rationale for 
the known inability of the skeleton expansion to capture such behaviour, and points to the intrinsic 
dangers of partial infinite-order summations that are based on PT in U. 

PACS numbers: 71.27.-|-a, 75.20.Hr, 71.30.+h. 



I. INTRODUCTION 

Since its inception nearly forty years ago, the An- 
derson impurity model (AIM)El has become one of the 
generic models of correlated electron physics, permitting 
at possibly the simplest level a detailed study of the com- 
petition between local interactions and electron itiner- 
ancy. The great majority of previous work, reviewed 
comprehensively in Ref. @, has focused on an impurity 
coupled to a conduction band whose width (D) repre- 
sents a large energy scale — large compared to the one- 
electron hybridization strength (Aq) that characterizes 
the host/impurity coupling. The opposite extceme of a 
narrow-band AIM has however been invokedaC^in rela- 
tion to the Hubbard model in large dimensionaj, where 
the correlated lattice problem maps onto an effective im- 
purity model coupled to a self-consistcntly determined 
bathQ. Here, in the vicinity of the Mott transition oc- 
curring in the particle- hole symmetric (i-fiUed) Hubbard 
model, an effective narrow-band AIM is generated self- 
consistently in the now widely accepted scenario for the 
transition^ that arose originally from the iterated pertur- 
bation theory (IPT) approach of Kotliar, Georges and 
coworkersQ'l. P 

Recently Hofstetter and Kehreinu have undertaken 
a numerical study of the symmetric narrow-band 
AIM, freed from the additional complications of self- 
consistency and itself defined by Aq ^ D. Thermody- 
namic and spectral properties of the model were found 
thereby to exhibit a range of, behaviour that is quite 
distinct from the conventionale wide-band AIM. Single- 
particle spectra for example, determined via finite-size 
Lanczos calculations, show a much richer structure than 
their wide-band counterpart. At low frequencies in par- 
ticular they exhibit an intricate competition between var- 
ious physical effects and their associated energy scales. 



reflecting ultimately the importance in the narrow-band 
AIM of an energy scale ojq ^ \/AoD that while large 
compared to the host bandwidth D is typically small 
compared to the on-site interaction U. 

Spectral properties of the narrow-band AIM are dis- 
cussed in the present paper, a complementary analytical 
study based in part upon 2nd order perturbation the- 
ory (2PT) in U. This is found to provide a qualitatively 
correct description on all energy scales, suggests a sim- 
ple physical interpretation of the most prominent spec- 
tral features, and produces rather good agre^ent with 
the Lanczos results of Hofstetter and Kehreiro. We also 
show, independent of the intrinsic limitations of 2PT, 
that spectral behaviour divides into two distinct regimes 
that are borne out by the calculations of Ref. ||: (a) 
weak-coupling for AqZ/D ^ 1 (with Z the quasi-particle 
weight), and (b) a strong coupling (Kondo) regime for 
AoZ/D < 1. 

For the weak to moderate interaction strengths reli- 
ably accf^ible by Lanczos calculations, Hofstetter and 
Kehreirjj'u have shown that pole contributions to the in- 
teraction self-energy Ei(a;), occurring on the scale ojq ^ 
VAo-D, are an integral facet of the narrow-band AIM; 
but that they cannot be captured to any order in a skele- 
ton expansion. Within 2PT, Sj(a;) is indeed found to be 
dominated entirely by poles on this order. However in 
strong coupling [/ ^ Aq, 2PT will certainly be qualita- 
tively deficient on the lowest energy scales \uj\ <^ D char- 
acteristic of the Kondo resonance. We show nonetheless 
that in strong coupling the pole contributions to the self- 
energy are again given precisely by the result arising from 
2PT. These poles are thus a ubiquitously dominant and 
essentially [/-independent characteristic of the narrow- 
band AIM. They are moreover continuously connected to 
the atomic limit of the model, which sheds light on the 
inability of a skeleton expansion to recover them, and 
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points to the intrinsic dangers of partial infinite-order 
summations based on PT in U. 

The Hamiltonian for the AIM is given in standard no- 
tation by 



+ 51 ("^'*kC+ Ck. + h.C.) (1) 

kcr 

where the first term describes electrons (of spin a =t, 4) 
in a metallic host band of dispersion ek- The following 
two terms refer to the impurity, with Si the impurity level 
and U the on-site Coulomb interaction. The final term 
describes the one-electron hybridization between the im- 
purity and host. Throughout this article, as in Ref. ^, we 
study the particle-hole symmetric AIM where £i ~ —U /2. 
In this case, for all interaction strengths, the Fermi level 
remains fixed at its non-interacting value and the impu- 
rity charge Ui = (n^i -I- ) = 1; single-particle spectra 
are thus symmetric with respect to the Fermi level, lu = 0. 



II. NON-INTERACTING LIMIT 

We first consider briefly the noninteracting problem 
Ei — ~U/2 = which, although rather trivial by it- 
self, generates two inherent energy scales that control the 
subsequent evolution of single-particle properties with 
interaction strength. For U = the (causal) single- 
particle impurity Green function G^{uj) = ReG°(w) — 
i sgn(tt')7r iD°(a;) is given bya 



A(c.)]- 
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with z = uj + iO+sgn(a;) and A{uj) — Aji{uj) — 
i sgn(cj)Ai(Lj) the hybridization function, A.{lu) — 

Ekl^.k|V(^-ek). 

We are interested in the narrow-band case, meaning 
D > Ao = Ai(w = 0) with D the host bandwidth. The 
physics of the model is naturally rather insensitive to the 
precise form of the hybridization. We thus take Ai{lu) to 
consist of a single flat band of intensity Aq , ranging from 
—D to +D: so that in total 



A(c^) = — log 
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where Ar(w) follows via Hilbert transformation. 

The noninteracting single-particle spectrum consists 
of two contributions: (i) a low-energy continuum for 
LU e [—D, +D], arising from the non-zero imaginary part 
of the hybridization; (ii) two poles, one lying above and 
the other one symmetrically below the continuum. The 
fraction of spectral weight residing in the poles depends 
strongly on the bandwidth of the host metal. In the usual 
wide-band model, -D Aq, most of the spectral weight is 



concentrated in the single-particle band, while the poles 
are exponentially weak and hence irrelevant in practice. 

The situation is however radically different in the 
present narrow-band model, defined hy D ^ Aq. Here 
the poles are the most prominent features of the spec- 
trum, suggesting that the narrow host behaves in effect 
as a single site or leveH. As this "host site" is coupled 
to the impurity, a bonding and an anti-bonding orbital 
form in analogy to the H2 molecule. These orbitals oc- 
cur at frequencies ±luq far outside the band, obtained by 
expanding Ar((j;) ~ (2/tt)AqD/lu in the denominator of 
(I) to yield 
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The pole- weight of each orbital isq — (2Ao — 7rD)/(4Ao — 
ttD) ~ ^ + O(i:>/Ao), showing that they carry almost 
the entire spectral density. 

The noninteracting single-particle spectrum is thus 
characterized by a weak continuum of weight O (D/ Aq), 
ranging over the lowest energy scale |a;| < D, and two 
molecular orbitals carrying most of the spectral inten- 
sity. The latter occur at the orbital bonding energy ojq 
which, although much larger than D, defines a second 
low-energy scale luq ^/ AqD <C Aq. Central ques- 
tions then arising include: how does this low-energy scale 
evolve with increasing interaction strength, over what 
range of U does it constitute the dominant low-energy 
physics of the problem, and when by contrast is the lat- 
ter dominated by the Kondo effect that one expects to 
prevail for sufficiently large interactions? It is these ques- 
tions we aim to address in the following sections. 



III. PERTURBATION THEORY 

The impurity Green function for the interacting sys- 
tem, Gi{u}), is related to G^{bj) by the Dyson equation 



(5) 



This simply defines the interaction self-energy Si(w) — 
Yi^{u}) — isgn(a;)I]J(a'), considered here via perturbation 
theory (PT) in U about the noninteracting limit. One 
does not doubt such an approach in principle for, al- 
though we are interested in the narrow-band regime, the 
impurity is coupled to a metallic host and one thus an- 
ticipates ubiquitously a Fermi liquid ground state th, 
is perturbatively connected to the noninteracting limiti 
PT is of course limited in practice, by the order to which 
it is taken. 

The 1st order (Hartree) contribution to E^, of = 
U/2, is cancelled precisely by the bare site-energy = 
—U/2. The lowest order process is thus of 2nd order 
in U, whereby a propagating tr-spin electron generates a 
particle-hole fluctuation of opposite spin on the impurity. 
Diagrammatically, it reads 
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with wavy lines representing the interaction U and left- 
going (right-going) solid lines representing a free particle 
(hole) propagator Gf. 



A. Particle-hole polarisation propagator 

The — (T-spin particle-hole bubble appearing in (|^) 
stands for the noninteracting polarisation propagator 

-t-oo 

ouiu:) J ^G'd^')GUco'-u;). (7) 

— oo 

Using the Hilbert transform relation for G^{lu), a short 
calculation gives 

+ \ui\ 

hm°n{uj) = 9 (lo) J duo' D^UJ') D1[lo' ~ Lo) 





+e{-Lo) J dc^'A"(cc;')A"(^^'-w) ■ (8) 

-\uj\ 

Since Im°n(a;) is symmetric with respect to the Fermi 
level w = 0, it is sufficient to discuss the first term 
(w > 0). ft describes particle-hole excitations obtained 
by transferring an electron from an occupied state be- 
low the Fermi level {uj' — o; < 0) to an unoccupied state 
above it (to' > 0). The narrow-band model is domi- 
nated by particle-hole excitations between the bonding 
and the anti-bonding orbital, involving an energy cost 
of 2u!o, which in ilm"n(ci;) translates to pole of weight 
~ 1/4 at w = 2uJo- 

In addition, Ini°n(ci;) contains two continua, associ- 
ated with particle- hole excitations to/from the single- 
particle band and therefore of small spectral weight. 
A first continuum arises for lu e [luq^luq + D] (net 
weight 0{D/Ao)), corresponding to electronic transfer 
between the single- particle band and one of the orbitals. 
A second appears for low frequencies |a;| < 2D (net 
weight O ([Z?/Ao]^)), arising from particle-hole excita- 
tions within the single-particle band. As D^{iu) is finite 
at the Fermi level lo = 0, the low- frequency behaviour of 
the latter contribution is 

ilm°n(c.) \u.\[D^{Qf (9) 

(which is used below to prove Fermi-liquid properties on 
the lowest energy-scale). 



B. Self-energy in 2PT 

The self-energy in 2nd order perturbation theory 
(2PT) is given by (^, which translates to: 

^d^) = J ^^"Tli^')G^{u;~co') (10) 

— oo 

In analogy to the above calculation for the polarisation 
propagator, we use the Hilbert transforms for "ll and G^ 
to obtain: 

mLj) = 0{Lo)U^ J dLu'Im°n{Lu')DUuJ-to') 





+e{-u;)U^ J dLu'lm°Il{Lu')D^iu;-u;') (11) 
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Again, due to particle-hole symmetry, we may restrict the 
following discussion to positive frequencies; and (w) 
follows from Sj(w) by Hilbert transformation. 

In the narrow-band limit, the 2PT self-energy is dom- 
inated by the following process: an extra cr-electron is 
placed in the empty anti-bonding orbital (energy cost ujq) 
and excites an electron of opposite spin from the bonding 
to the anti-bonding orbital (energy cost 2ujo)- In ^K^); 
this generates a pole of weight ~ ttU'^/S at cj = Sujq, that 
gives by far the dominant contribution. Much weaker 
contributions to Sj(a;) arise in addition from the ex- 
tra electron being placed in the single-particle band in- 
stead of the orbital, or involve one of the particle-hole 
continua discussed above. This yields self-energy con- 
tinua for UJ e [2ujn, 2cjo + D] (net weight O (J/^D/Aq)), 
Lo e [ujo,LUo + 2D] (net weight O {U'^[D/ Aq]'^)), and 
\uj\ < 3D with net weight O (U^jD/Aof). The latter re- 
sults from a convolution, as in (fTl), of the low- frequency 
continuum of Im°H(a;) oc oj (see eq. (||)) and the single- 
particle band. It leads to Sj(w) oc close to the Fermi 
level to =jO, corresponding to conventional Fermi liquid 
behavioum. 



C. Spectral evolution in the narrow-band regime. 

Numerical results for the impurity spectrum Di{uj) = 
— ^agn{uj)lmGi{ijj) within 2PT will be given in the fol- 
lowing section. Here we develop an approximate, but 
rather accurate analysis that shows in simple physi- 
cal terms how the spectrum evolves with interaction 
strength, particularly on the low-energy scales of primary 
interest. 

As discussed above the self-energy is dominated by 
poles at UJ — ±30^0, its continua being negligibly small 
in the narrow-band regime and vanishing as I? — > 0. We 
thus retain solely the pole contributions, i.e. 
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with a corresponding quasi-particle weight Z 
{dY.f{uj)/duj)uj=oY^ given by 



Z 



U 

6^0 



H -1 



(12) 
[1- 

(13) 



The evolution oi Z = Z{U) is thus controlled simply by 
the ratio U /uja, and for U uji^ ^ \/Ao-D in particular, 



U 



72 AoL* 



(14) 



With Si(aj) given from (p^, Gi{ijj) follows directly 
from (H) . The resultant spectrum then consists of a con- 
tinuum for I a; I < D, and two sets of poles at |w| = lo+ 
and oj- ( "set" meaning a pole at positive frequency and 
its counterpart occurring symmetrically below the Fermi 
level Ld — {)). The latter follow from the single-particle 
pole equation w — Ar(w) — Sf^(w) = 0; insertion of (p^ ) 
into which, together with the expansion AR(a;) ~ '^^^ 
— which is valid for 3> D — yields a biquadratic with 
solutions 
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(15) 



The set of low- frequency poles at cj- < wqi with pole 
weights (J-, have their origin in the molecular orbitals 
of the non-interacting problem, to which they reduce for 
U/ujQ where — > ujq and g_ 1/2. In addition, 
for any [/ > 0, a second set of poles arise at a;+ > 3aJo- 
In weak coupling [//wq ^ li the latter have negligible 
pole-weight, g+ ~ {U/ojo)'^. With increasing ?7 however, 
they drain weight from the low-frequency orbital rem- 
nants and shift to progressively higher frequencies until, 
for U ^ ujQ ^ ^/ AqD, they correspond to the Hubbard 
satellites and overwhelmingly dominate the spectrum: 
UJ+ - U/2 with pole-weight q+ - 1/2 - {[ivo/U]'^)- 

The above pole analysis is exact in the strict limit of 
vanishing bandwidth D 0, for fixed ujq ~ \/ AqD. 
Here, the host behaves precisely as a single orbital at the 
Fermi level w = 0, and the problem reduces to the sim- 
ple two-site model discussed by Langal3 and employed 
by Hofstetter and Kehreincl to rationalize their Lanczos- 
determined spectra for the narrow-band AIM. In this 
limit the hybridization Ai(uj) is a ^-function, and its 
weight 
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relates the bonding energy loq (eq. (^) to the hopping 
between host and impurity, loq — Vio. As pointed out 



by Lang , 2PT is in fact exact in this case, the 4-pole 
spectrum being given precisely by eq. (|l^); for U ^ uiq 
in particular, the low-energy pole lj- reduces to 
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where J = AV^/U is the antiferromagnetic exchange cou- 
pling between the impurity and the host orbital. 

As a "zero-bandwidth" approximation, the two-site 
limit just discussed is however rather limited as a means 
to interpret the narrow-band AIM. It is certainly accu- 
rate for the high-energy spectral poles uj+ (eq. ([l5|)) that 
evolve for U ^ utQto the Hubbard satellites at 0-1+ ^ U/2. 
But it naturally fails both on the lowest energy scale of 
the bandwidth D and, relatedly, in capturing the strong 
coupling behaviour of the low-energy spectral poles at 
= LU-: the latter, by definition, can only occur out- 
side the band (u;_ > D), whence eq. ( p7| ) suggests a 
breakdown of the two-site analogy for J ^ O (D) or 
U^OiAo). 

It is however straightforward to obtain, for all in- 
teraction strengths, a uniform description of the spec- 
trum at low energies; where by low energy we mean 
u! < ujQ ^ \/ AqD, thus encompassing both the lower- 
energy poles w_ and the continuum for |w| < D. For 
provided only that |ti;| <C 3a;o, the self-energy given 
from ( p^ ) reduces to its leading low-frequency expansion 
Sf (w) ^ — (1/Z — l)ti^- Hence, for |a;| <C 3u;o, the impu- 
rity Green function Gi{uj) is given from (^), using (||) for 
the hybridization, by 



AoG,(c^) = 
— — log 

Aoii TT 



(18) 
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where Co — lo/D. This equation may be used to under- 
stand the low-energy spectrum for all U . It depends upon 
the single parameter AqZ/ D, where the {/-dependence of 
the quasi-particle weight Z is given explicitly by eq. ( |T^ ) 
for 2PT. And the characteristic spectra divide into two 
distinct regimes, according to whether AqZ/D ^ 1 or 
^ 1, as now considered. 



1. AoZ/D > 1. 

From eq. (|l|), AqZ/D > 1 corresponds to CZ/Aq < 
^^J^Jtx ~ 5. Here the spectral poles, at = occur 

far outside the band. Specifically, expanding log 



2/lo for PI > 1, eq. 



yields = ./l^ » 1, i.e. 



(19a) 



with pole-weight 
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= 2 



(19b) 



Thi s is essentially the two-site limit result0 for — 
eq. (19a) with Z from (|l^) differs negligibly from eq. ( [l5| ) 
— encompassing both weak coupling U /loq where 
Z — > 1 and uj^ — > loq] as well as J7 ^ ~ VAqD 
(subject to U/Ao ^ 5), where uj^ reduces to eq. (17). 
In addition of course, eq. ( p^ ) yields correctly a spectral 
continuum for |u;| = \uj\/D < 1. But since ^^I'^l ^ ^ 
for all \u}\ < 1, the continuum reduces to that character- 
istic of the non-interacting limit, U = (where 
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U! m 



eq. (|18|) is likewise negligible for the narrow-band model). 

Hence, for AqZ > D (i.e. U/Aq < 5), the spectral 
'action' takes place exclusively in the polej;,^ (and cl;_|_) 



and is well captured by the two-site limi while the 
low-energy continuum is essentially unaffected by elec- 
tron interactions. 



2. AoZ/D < 1. 

The above picture is marked contrast to that arising 
from AqZ <C D (i.e. U/Aq ^ 5), where the uj- poles are 
exponentially close to the band edge, being given from 



(N)by 



2 exp 
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The corresponding pole-weight is however exponentially 
small. 



The essential conclusions above are nonetheless largely 
independent of 2PT. Eq. ([l^ ) amounts to a low-frequency 
expansion of the self-energy (embodied in the first term) 
that is quite general, independent of the {/-dependence of 
Z, and valid for ^j'l'^l ^ ^- This by itself is sufficient to 
infer the two distinct spectral regimes, AqZ 3> D where 
the spectrum contains low-energy poles at lo^ (given in 
terms of ujq and Z by eqs. ([l^)), and with a continuum 
for < D that is essentially unrenormalized by in- 
teractions; and AqZ <C D where solely the continuum 
survives and eq. ( pT| ) provides a reasonable description 
of the Kondo resonance for uj/Aa2 <c 1. Further, since 
2PT is exact for the two-site limitlij, we anticipate that it 
should prowde a sound description of the observed Lanc- 
zos spectraB at least for AqZ/Z? > 1, as considered in the 
following section. 

Before proceeding we comment, briefly on the sug- 
gestion by Hofstetter and KehreinH that the low-energy 
spectrum (including the poles) in the narrow-band limit 
D ^ Q, should be described by a scaling function 
Di{uj) = f{uj — Lo/D) where / is independent of 
D. Eq. ([ij) supports this, since from the 2PT result 
(|l|) AqZ/D ^ [Ap /Uf is independent of D provided 
U ':$> u)Q ^ \/AqD (which encompasses both the 
AqZ/ D ^ 1 and ^ 1 regimes). This however is distinct 
from the one-parameter universal scaling that is charac- 
teristic of the strong coupling Kondo regimeH. The latter 
arises only for AqZ <C Z?, where the low-energy spectrum 
Di(uj) = f{uj/ AqZ) but is otherwise independent of any 
of the bare parameters D, Aq or U. 
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so the low-energy poles in effect are irrelevant. 

All the action now takes place in the continuum: from 
(p^), and for all < 1 save exponentially close to the 
band edges, the spectrum for AqZ <^ D is given by 



7rAoA(w) = 



AqZ 



(21) 



This of course is the Kondo resonance, scaling univer- 
sally in terms of the single Kondo scale uik = AqZ but 
with no explicit dependence upon the bare parameters 
D, An or U. Its Lorentzian shape, arising also (see eg. 
Ref. y) in microscopic Fermi liquid theory and slave bo- 
son approaches, is but a caricature of the true Kondo 
resonance — which is known for example to contain long 
(Doniach-Sunjic) tailsM'll3 that are not captured by the 
simple Lorentzian form. Moreover the [/-dependence of 
the quasi-particle weight Z, given explicitly by eq. ( p^ ) 
within 2PT, fails to capture the exact asymptotic be- 
haviour Z cx exp(— 7rC//(8Ao)) appropriate to the strong 
coupling Kondo regime U /Aq 3> 1, producing instead an 
algebraic decay with U for all U ^ ujq ^ ^/ AqD as in 



D. Spectra: numerical results 

We turn now to the numerically determined 2PT spec- 
tra. Fig. shows the low- frequency behaviour of Di{uj) 
vs. u! — oj/D, for a bandwidth of I? = 10^''Ao and 
interaction strength U = 0.2Ao. For these parame- 
ters AqZ/D ~ 540, safely in the weak coupling regime 
AqZ/D ^1 — albeit 'strong' coupling in the sense 
that U/luq ~ 25, as manifest in a quasi-particle weight 
Z ~ 0.054 <C 1 that is far removed from the non- 
interacting limit where Z = 1. The spectrum is shown 
for \uj\/D < 200, and we remind the reader that almost 
all (~ 95%) of the spectral weight resides in the Hubbard 
satellites at frequencies as large as U/2 = 10"^ I?, which 
are omitted from the figure. 

On the lowest energy scale \uj\ < D we see the ex- 
pected continuum of net weight O (D/ Aq) which, in ac- 
cordance with the discussion above, is indeed essentially 
unrenormalized from the non-interacting limit. At the 
band edges lu — ±D, the spectrum vanishes of necessity 
since the hybridization A^{tu) is singular there. Due to 
the lowest frequency self-energy continua (§ [II B above), 
the central continuum acquires side bands ranging from 
D < \uj\ < 3D, whose weight remains tiny over the whole 
range of interactions. 
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On the next higher energy scale, proportional to the 
antiferromagnetic coupling J, we find the 'orbital rem- 
nants' at (eq. 0) lj- = 3J/2 = U/^^D / {ttU) ~ 1%D. 
These poles are by far the most prominent features of 
the low-energy spectrum, carrying almost all the spec- 
tral weigh t in this sector, although their absolute weight 
(eq. ( |l9b| )) = Z/2 ~ 0.027 is small since the domi- 
nant spectral intensity is carried by the high-energy Hub- 
bard satellites. A third low-energy scale is the molecular 
bonding energy uoq ^ \J AqI? emerging from the non- 
interacting limit. It enters th e sing le-particle spectrum 
via the self-energy continua (§ IIIB) and yields two sets 
of weak bands at |w| ~ — SOD and |w| ~ Iloq ~ 1601?, 
albeit with net spectral weights that render them insignif- 
icant in practice. 

A comparison of Fig. ^ with the spectrum of an 11 + 1- 
site Anderson star, calculated numerically for the same 
parameters by Hofstetter and Kehreiio with the Lanczos 
method, shows that 2PT is able to reproduce essentially 
all relevant energy scales. Good agreement is found in 
the lowest energy scale D, given by the Fermi liquid con- 
tinuum, and on the high-energy scale t//2 for the Hub- 
bard satellites, although 2PT remains somewhat incom- 
plete on the antiferromagnetic scale J. Here, instead of 
the single pole predicted by 2PT at \uj\ ~ 19Z3, 

the Lanczos calculations show two sets of (5-peaks, the 
first occurring for |a;| ~ 18 — 2113 and the second for 
I a; I ~ 25 — 26Z3. Although it is not clear how these sets 
will evolve with a larger number of sites in the Lanczos 
calculations, a structure with at least two features will 
most likely be preserved in this region. This points to 
the natural limitations of 2PT, and the need for more 
sophisticated many-body techniques. 

For obvious numerical reasons, the Lanczos method 
cannot detect any features as weak as the 2PT con- 
tinua at loq ~ SOD and IijJq. Nevertheless we believe 
these continua to be robust. They cannot be cancelled 
by higher order self-energy diagrams since, for any fre- 
quency, > for all self-energy contributions. 

Having discussed the weak coupling regime /^.^Z/ D ^ 
1, we now turn to strong coupling, I^qZ/D <^ 1, and 
the approach to it. Fig. shows the low-energy Diiuj) 
vs. Cj = uj / D, again for a bandwidth of -D = 10~^Ao; 
and for interaction strengths U/ttAq = 2,4 and 6, corre- 
sponding respectively to AqZ/D ~ 0.58, 0.15 and 0.06. 
In agreement with the analysis of § [II C above, and ex- 
cluding the Hubbard satellites at ztU /2, the spectra in 
all cases are in practice confined exclusively to the low- 
frequency continuum |w| < D, whose persistent Fermi 
liquid character is manifest in the pinning of the spectra 
at the Fermi level (uj = 0) to its non-interacting value 
of I/ttAq. For the two larger interaction strengths in 
particular, the Lorentzian character of the 2PT Kondo 
resonance arising for AqZ/D <C 1 is seen clearly. 

The results above concm; with the Lanczos calculations 
of Hofstetter and KchrcintI, performed at a moderate in- 
teraction strength U/Aq = 4. No low-energy poles arose, 
and the spectrum was found to be governed by a Fermi 



liquid continuum with a broad resonance centered on the 
Fermi level and two sharper features at the band edges. 
As shown by the top curve of Fig. ^, 2PT reproduces 
this behaviour remarkably well, albeit that it sets in at a 
somewhat larger interaction strength of U/ Aq ~ 6. 

For numerical reasons, the strong coupling Kondo 
regime AqZ/D <C 1 naturally cannot be captured ad- 
equately by the Lanczos calculations^, and the deficien- 
cies of 2P T in this regime have also been highlighted 
It 



[HC 



appears however that 2PT provides a 
rather reasonable description of the Lanczos results for 
the narrow-band AIM, and is valid up to interaction 
strengths U/ttAq on the order of 2 or so, a regime of 
validity that corresponds closely to that found also with 
2PT for the normal wide-band (D = oo) model (see eg. 
Ref. |). 



IV. SELF-ENERGY, AND SKELETON 
EXPANSION. 

Within 2PT the structure of the self-energy for the 
narrow-band model (Aq S> D) is rather simple, being 
dominated by poles at w = ±3wo: each of weight 
(see eq. (|T2|)). In view of the natural limitations of 2PT 
however, one might ask what may be inferred more gener- 
ally about the behaviour of Ei(tj), and in particular pole 
contributions thereto? It is aspects of this issue that we 
now consider. 

The self-energy is effectively defined via the Dyson 
equation (H) whence, for frequencies lo where the single- 
particle spectrum is purely real, Si(w) is given by 



Ei(cj) = UJ - AR(tj) - 



1 



(22) 



Here X{uj) = KeGi{uj) is thus defined for convenience, 
and is given by the Hilbert transform 



X{uj) = 2lu du 



(23) 



(where particle-hole symmetry has been employed, and 
a principal value is implicit). As is well knownd (see also 
Ref. it follows from eq. ( p^ ) that the zeros of X(uj) 
correspond to the poles in Si(aj). For weak to moder- 
ate intecaction strengths {U/Aq < 4), Hofstetter and 
Kehreim find that the Lanczos determined X{uj) = 
for UJ ^ O {ujq) and hence that 'Si{uj) contains poles at 
frequencies on the order of ±\/ AqD. This concurs with 
2PT, which we have argued (§ IH) to be valid in such 
a weak coupling domain, where AqZ/D > 1. But what 
of the strong coupling regime U ^ Aq (3> D), where 
AqZ/D <C 1 and the Kondo effect is well developed? The 
limitations of the Lanczos technique in this domain have 
been alluded to above, but the general structure of the 
Lanczos-determined spectrum Di{u>) is in fact sufficient 
to determine the poles of Ej(a;), as now shown. 
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In [Strong coupling the single-particle spectrum is 
founda to separate into high- and low-energy parts, 
= Df{Lj) + Df{uj); we focus below on > (by 
particle- hole symmetry). The high-energy contribution 
is Df-{uj) — q+S{uj — 0J+)- It corresponds to the Hub- 
bard satellite, with position a;-|_ ^ U/2 and pole-weight 
5+ ~ 1/2. By contrast the low-energy continuum Dfloj) 
is non-zero only on the band scale |a;| < D, is pinned at 
the Fermi level where 7rAo-Di(a' 0) = 1 ubiquitously, 
and contains a well developed Kondo resonance on the 
scale UJ ^ AqZ <^ D. From eq. ( p3| ) it follows that 
X(lu) = Xii{ll!) + Xi^{lu), given asymptotically by 



u;«C//2 2g+ 

Xh(w) - 772 



Xl{uj) 



(24a) 
(24b) 



where 6 = dcjDfluj) is the integrated weight of the 
low-energy continuum. From eq. (|^ the zeros of X{uj), 
and hence the poles of Y,i(uj), thus occur at cj = ±LUp 
given by 



UJp = U! 




(25) 



(which is self-consistent provided D <C ^ U/2 as in- 
deed found, see below); the corresponding pole- weight in 
i;,(w), given from eq. (||) by Q = TT/\dX{ui)/dLj\^^^^, 
is 



Q = 



(26a) 



Since Di{uj) is normalized to unity however, the inte- 
grated low-energy spectral weight 5=1 — 25+ . In strong 
couphng moreover, where 2\ei\ = U » Aq (S> D), the 
leading asymptotic behaviour of the high-energy spectral 
pole-weight (7+ is in turn readily deduced by taking the 
limit of constant Vi\^ and D ^ 0, i.e. by considering fixed 
loq ^ VAe-P a-nd — > 0, which is just the two-site limit 
for whichEJ q+ ~ 1/2 — 18a;§/C/^ (amounting in effect to 
perturbation theory in Vi\^/U). Hence 5 ^ SGwq/L/^, and 
from eq. ^ 



UJp 



3wi 



• 



(26b) 



Eq. ( p6b| ) (with Q from ( p6a| )) is the result in strong 
coupling U S> Aq for the poles in Y,i{uj). The arguments 
for it do not depend upon perturbation theory in U about 
the non-interacting limit, but the result itself is precisely 
that arising from simple 2PT, eq. (|l2[). This is salu- 
tory, for although 2PT extrapolated to strong coupling 
produces an inadequate caricature of the Kondo reso- 
nance on scales w ~ O (Aq.^) <C D (see e.g. § HIC), the 
above arguments suggest that on scales of order wo ~ 
V Aq-D it becomes asymptotically exact in strong cou- 
pling. In addition, noting that ujq is [/-independent and 



that eq. (26b) holds both perturbatively for U ^ O (Aq) 
and in strong coupling [/ ^ Aq, it is clear that any de- 
pendence of the poles in Si(w) upon U is at best weak 
across the entire range of interaction strengths. 



A. Skeleton expansion 



There is a second feature of the above analysis to which 
we would draw attention, namely that the poles in (uj) 
connect continuously to that characteristic of the atomic 
limit, Aq = (i.e. Viu = 0). As Ao the two poles 
in Ei(w), each of weight Q = ttC/^/B and occurring at 
UJ = ±3ujo with ujQ ~ ^/ AqD, coalesce to a single pole at 
UJ = with weight 7rJ7^/4, thus producing 



^t{uj) 



1 



A UJ + iO+sgn(w) 



(27) 



This result, which corresponds to 2PT, is of course well 
knowio to be exact in the atomic limit for all f/ > 0, a 
feature that is particular to the particlc-holc symmetric 
case considered. 

Hofstetter and Kehreini have considered the behaviour 
of Tji{uj) for the narrow-band AIM from the perspective of 
the skeleton expansion, where the exact propagator Gi (uj) 
is inserted into every skeleton diagram contributing to 
Si(w). They point out that the pinning of the single- 
particle spectrum at the Fermi level (for all C/ > and 
any Aq > 0) attests to the convergence of the skeleton ex- 
pansion on the lowest energy scales characteristic of tlie 
Fermi liquid continuum. Conversely, following KehrcinQ, 
they show that the |a;| ^ O {ujq) poles in 'Si{uj) cannot be 
explained in any order of the skeleton expansion, which 
fails to converge for such frequencies. Given the above 
connection between the poles in Si(w) and that endemic 
to the atomic limit, the latter behaviour is entirely natu- 
ral, since it is known that the skeleton expansion fails to 
converge for the atomic limit (and miare generally, but in 
essence equivalently, for an insulatora). 

In view of the above it may be instructive to con- 
sider the failure of the skeleton expansion, taking the 
atomic limit as a paradigm. The essence of the skeleton 
expansion to any finite order is that it involves a par- 
tial infinite-order summation of diagrams obtained from 
PT in U (and themselves expressed in terms of the non- 
interacting propagator G^(tj)). As such however, it may 
fail to include higher-order diagrams that act in large 
part to cancel those that are included, and this is where 
the dangers arise. 

Consider for example the following diagram, which is 
4th order in PT but is included in the class of 2nd order 
skeletons: 
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(28) 



It represents a process creating two particle-hole excita- 
tions of spin —a on the impurity site. If the first such 
pair has not hopped off the impurity when the second 
pair is created — as is inexorable in the atomic limit — 
then such processes are formally forbidden by the Pauli 
principle. They are nonetheless properly cancelled by ex- 
change diagrams arising from the same order in PT, for 
example: 




(29) 



However the exchange diagram (^9|), while arising to the 
same order in PT as (p8|), belongs to the class of 4th 
order skeletons. In this sense it is of higher order, and 
would not therefore be included if the skeleton expansion 
was truncated at 2nd order. 

Generalizing this reasoning illustrates that the skeleton 
expansion, if truncated after a certain class, will include 
specific series of diagrams without necessarily accounting 
for the corresponding exchange terms, since the latter be- 
long in general to a higher skeleton class. Hence, even for 
systems where PT in U is known to converge order-by- 
order, the order-by-order convergence in skeletons may 
not be taken for granted. 

The atomic limit, whose relevance to the narrow-band 
AIM has been pointed out above, provides a direct ex- 
ample of the latter point. Here Si(w) is given exactly by 
2PT (eq. (|2^)), i.e. by I]i(w) = ^G°{uj) (where, triv- 
ially, the unperturbed G^{uj) ~ l/{uj + iO+sgn(a;)). The 
problem is convergent order-by-order in straight PT, all 
diagrams contributing to Si(w) in any given order n > 2 
of PT summing precisely to zero. Consider by contrast 

(2) 

the second-order skeleton expansion, denoted Sg^ei and 
obtained from eq. (^) by replacing G° with the exact Gi , 

togive01]^^ii(c^) = ^u/[u^-{W/2)^] -u/9. This 
fails entireiy to capture the ^ 1/lu behaviour of the ex- 
act Ei((jj)ll3 (save trivially for |ti>| 3> 3U/2), and that it 
does so for the essential reasons outlined above is directly 
evident by recasting Sg^gj as a functional of G^, viz: 



.(2) 
■^skcl 



1 



[fGO(^)]- 



(30) 



Finally, we note that the behaviour described in this 
section is not specific to the narrow-band AIM, but arises 



also in the infinite-dimensional Hubbard model at half 
filling (in which context study of tiia-Jiarrow-band AIM 
was first motivatedB). Within IPTETH, as U approaches 
the critical Uc for the Mott transition from the metallic 
phase U < Uc, the self-energy acquires poles (or strictly, 
sharp resonances) at finite a; in a preformed gap; and as 
U Uc— the pole positions approach the Fermi level 
uj = (less rapidly than the central Fermi liquid con- 
tinuum vanishes), and with weights that remain finite. 
Hence as C/ — > Uc— the low-frequency behaviour char- 
acteristic of an insulator — Si(tj) ~ A/{u! + iO+sgn((jj)) 
which amounts in essence to that for the atomic limit — 
is smoothly recovered via IPT. The poles in Si(w) are 
thus an integral facet of the Mott transition and, from 
the discussion above, we believe the success of IBX in 
capturing it is intimately connected to its abilityEl^B to 
recover correctly the atomic limit. By contrast, since 
the poles in Si(w) cannot be captured via the skeleton 
expansionQ, approaches based at heart upon|-tljie latter 
— such as self-consistent perturbation theoryt^ — sim- 
ply fail to uncover the transitionEj. 



V. CONCLUSIONS 

We have studied in this article the symmetric AIM 
with a narrow host band, £> <C Aq. Simple though it is, 
second order perturbation theory in U is found to give 
a rich and relatively complete account of the underlying 
single-particle dynamics. In particjular, and in agreement 
with recent Lanczos calculationso, it leads naturally to 
two distinct regimes of spectral behaviour according to 
whether AqZ/D 1 or ^ 1, and whose essential char- 
acteristics we have argued to be largely independent of 
the details of 2PT itself. 

We have also shown that 2PT is remarkably robust 
on energy scales of order ~ \/ AqD, that reflect the 
underlying molecular orbitals characteristic of the non- 
interacting limit. Here we have argued that the 2PT re- 
sult for the self-energy is correct both perturbatively in 
weak coupling (by construction), and in the strong cou- 
pling limit U 3> Aq. Such poles in Si(u;) — which ca 
not be captured to any order in a skeleton expansion! 
— thus appear to be an essentially {/-independent char- 
acteristic of the narrow-band AIM, and that they are 
captured by 2PT is in turn closely related to its ability 
to recover correctly the atomic limit of the model. 

The latter comments also illustrate the difficulties 
that confront any theory which seeks to describe single- 
particle dynamics of the narrow-band AIM in strong 
coupling, and on all energy scales. At low energies 
UJ ~ llik ^ D characteristic of the Kondo resonance, low- 
order PT in U will naturally not suffice and an intrinsi- 
cally non-perturbative approach will be needed. And yet 
for energy scales on the order of ^/ AqD ^ Ag the essen- 
tial result of 2PT must be recovered — which as we have 
discussed is a delicate matter, since 'obvious' approaches 
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based on partial infinite-order summation of PT in U 
are liable to qualitative failure here. While accessible to 
numerical techniques such as the numerical renormaliza- 



retical approach that can encompass these twin dictates. 
In a subsequent paper we will however Apw that a re- 



the problem. 
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FIG. 1. Low-frequency 2PT impurity spectrum, for band- 
width D = 10"*Ao and U = 0.2Ao, where AqZ/D > 1. 



Note the logarithmic scale for Di{uj) 
The Hubbard satellites at |w| ~ U/2 -- 
Pull discussion in text. 



employed for clarity, 
are not shown. 




FIG. 2. Low-frequency 2PT impurity spectrum, 7rAo/)i(w) 
vs. oj/D, for bandwidth D = 10"" Ao and U/Ao = 2n (top), 
47r (middle) and Gtt (bottom), corresponding respectively to 
AoZ/D ~ 0.58, 0.15 and 0.06. The emergence of the Kondo 
resonance for AqZ/D <C 1 is seen clearly. Full discussion in 
text. 



10 



